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1.0  INTRODUCTION 

1.1  General 


The  knowledge  of  the  hydrologic  conditions  which  exist  during  low 
flow  periods  can  be  of  primary  importance  in  undertaking  a  variety 
of  water  resources  investigations.  When  analysing  water  quality 
conditions,  the  low  flow  characteristics  within  the  watercourse  are 
a  major  concern  to  both  the  user  and  the  Ontario  Ministry  of  the 
Environment.  Some  uses  of  low  flow  information  include  the 
fol lowing: 

i)  The  analysis  of  municipal  and  industrial  effluent  discharg- 
es to  streams 
ii)  Instream  pollutant  analyses  (point  and  non-point  sources 
Including  mixing  zones) 
iii)  Reservoir  design  (low  flow  augmentation) 
iv)  Environmental  approvals 
v)  Feasibility  of  small  hydro  developments 
vi)  Water  supply  and  evaluation  for  water  taking  permits 
vii)  Base  flow/groundwater  recharge  and/or  contamination  analy- 
sis 
viii)  Stream  fisheries  assessments 

ix)  Analyse  effects  of  changes  in  watershed  on  low  flows  (e.g. 

deforestation,  urbanization) 
x)  Agricultural 
x1)  Other 

The  identification  of  low  flow  characteristics  within  a  watercourse 
is  most  easily  accomplished  using  continuous  hydrometric  data 
recorded  for  the  stream.  Available  data  were  previously  analysed 
and  published  as  a  series  of  "Low  Flow  Characteristics"  maps  which 
were  published  several  years  ago  (more  than  10  years  In  most  cases) 
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by  the  Ministry  of  the  Environment  (Ministry  of  the  Environment, 
1973-1978).  Additional  data  (both  temporally  and  spatially)  are 
now  available  and  it  was,  therefore,  considered  appropriate  to 
update  the  available  data  base  describing  low  flow  characteristics 
across  the  Province  of  Ontario. 

1.2    Study  Objectives  and  Scope 

The  overall  objective  of  the  proposed  investigation  was  to  carry 
out  an  analysis  describing  the  low  flow  characteristics  at  suitable 
Water  Survey  of  Canada  streamflow  locations  in  Ontario.  The  foll- 
owing points  summarize  the  scope  of  the  investigations: 

1.  Update  statistical  low  flow  analyses  for  each  of  the  follow- 
ing administrative  regions: 

-  Central 

-  Southeastern 

-  Southwestern  and  West  Central 

-  Northeastern 

-  Northwestern 

2.  Undertake  data  base  screening  analyses  to  identify 
constraints  on  the  data  base 

3.  Produce  extreme  value  analyses  for  suitable  stream  gauging 
locations  (1,  3,  7,  15,  30-day  durations)  with  greater  than 
10  years  of  data 

4.  Produce  annual  flow  duration  analyses,  using  daily  data  for 
suitable  stream  gauging  locations 

5.  Produce  monthly  low  flow  analyses 
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2.0  METHODOLOGY 

2.1  General 


The  available  streamflow  data  base  was  obtained  for  all  relevant 
Water  Survey  of  Canada  monitoring  locations  (see  discussion  in  Sec- 
tion 2.2).  Data  analyses  and  screening  were  then  undertaken  in 
order  to  assess  the  usefulness  of  the  data  at  each  location  prior 
to  undertaking  statistical  analyses  (see  Section  2.3).  Extreme 
values  for  selected  low  flow  durations  were  then  calculated  for 
various  recurrence  intervals  for  both  annual  and  monthly  series 
using  the  techniques  discussed  in  Section  2.4.  Flow  duration  anal- 
yses were  also  undertaken  utilizing  the  daily  data  base  on  both  an 
annual  and  a  monthly  basis  (see  Sections  2.5  and  2.6).  In  most 
cases  the  extreme  value  analysis  results  should  be  used  for  low 
flow  characteristic  analysis.  However,  under  certain  conditions 
flow  duration  results  of  heavily  regulated  stations  may  give  more 
conservative  results.  Computer  drawn  maps  depicting  low  flow 
characteristics  were  then  produced  in  order  to  summarize  selected 
low  flow  statistics  on  a  regional  basis  (see  Section  2.7). 


2.2    Data  Base 


The  entire  Water  Survey  of  Canada  daily  streamflow  data  base  for 
the  Province  of  Ontario  for  the  period  of  record  to  the  end  of  1986 
was  obtained  in  computer  tape  format. 

Stations  with  a  minimum  of  ten  years  of  record  were  considered  to 
have  an  appropriate  record  length  for  the  purpose  of  this  investi- 
gation. Stations  Inactive  prior  to  1981  were  not  considered  In  the 
analysis. 

Several  data  management  programs  developed  by  Gumming  Cockburn 
Limited  were  then  used  for  extraction  and  analysis  of  data  from  the 
computer  archives.  First  a  program  was  written  and  executed  to 
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check  and  screen  the  available  data  base  at  each  station  for  miss- 
ing data.  A  second  program  was  then  run  to  utilize  this  analysis 
to  compute  running  averages  for  various  durations  and  to  pick  the 
corresponding  minimum  annual  low  flow. 

A  third  program  sorted  and  analysed  the  available  data  to  determine 
low  flow  values  on  a  monthly  basis.  Average  flows  were  determined 
and  extracted  for  the  annual  consecutive  low  1,  3,  7,  15  and  30  day 
durations  and  are  available  as  part  of  the  background  files. 
Another  program  was  used  to  compute  flow  duration  curves  as  dis- 
cussed in  Section  2.5. 

2.3   Data  Analysis  and  Screening 

The  data  were  manually  screened  to  remove  station  records  which 
contained  a  significant  number  of  zero  low  flow  occurrences  for  all 
durations.  These  stations  are  listed  in  Table  1  and  the  screening 
results  are  discussed  in  Section  3.2. 

The  initial  version  of  the  Low  Flow  Frequency  Analysis  program  was 
obtained  from  Environment  Canada  for  use  in  the  investigation 
(Pilon  and  Jackson,  1988).  This  program  contains  a  set  of  data 
analysis  and  checking  modules  suitable  for  data  screening  and  anal- 
ysis (Pilon  et  al ,  1985).  At  present  there  is  some  question  con- 
cerning the  reliability  of  the  methodology  for  checking  Independ- 
ence for  low  flow  data  series  (personal  communication,  P.  Pilon, 
Water  Resources  Branch,  Environment  Canada). 

The  data  screening  process  was  undertaken  to  test  the  assumption 
that  the  data  are  reliable  measurements  and  are  mutually  exclusive 
and  represent  independent  random  events  which  are  free  from  trend. 
The  non-parametric  tests  applied  include  the  Spearman  Rank  Order 
Correlation  Coefficient  for  Independence;  the  Spearman  Rank  Order 
Correlation  Coefficient  for  Trend;  and  a  general  randomness  test. 
The  statistical  tests  were  applied  at  the  1%  and  57.  levels  of  sig- 
nificance. (See  Appendix  A. 3  for  more  details  on  the  statistical 
tests.) 
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It  must  also  be  expected  that  a  certain  number  of  stations  which 
are  tested  will  randomly  "fail"  any  given  statistical  test.  There- 
fore, the  binomial  distribution  was  used  to  calculate  the  number  of 
expected  random  failures  over  the  entire  data  base  for  each  region 
in  order  to  assess  the  suitability  of  the  low  flow  data  on  a 
regional  basis.  A  tabulation  of  theoretical  and  actual  "failure" 
results  is  summarized  and  discussed  in  Appendix  A. 3. 4. 

The  effect  of  regulation  was  examined  by  further  subdivision  of  the 
data  base  according  to  regulated  and  non-regulated  stations.  These 
results  are  tabulated  in  Appendices  B  to  F  and  discussed  in  Section 
3.2. 

2.4    Extreme  Value  Analysis 

An  extreme  value  analysis  was  undertaken  for  each  of  the  1,  3,  7, 
15  and  30-day  durations  for  each  of  the  stations. 

For  example,  the  minimum  annual  consecutive  7  day  average  low  flow 
was  determined  for  each  year  for  the  station  record  being  analysed 
(see  Section  2.2).  The  corresponding  set  of  consecutive  7  day 
average  low  flows,  therefore,  represents  an  extreme  value  series  to 
which  a  theoretical  extreme  value  distribution  can  be  fit  for  the 
purpose  of  determining  the  low  flows  corresponding  to  various 
recurrence  intervals.  For  example,  the  7Q20  low  flow  is  often  of 
interest  for  water  quality  investigation  (i.e.  the  consecutive  7 
day  average  low  flow  with  an  average  recurrence  interval  of  once  in 
20  years). 

Previous  Investigations  have  found  that  the  Gumbel  III  distribution 
has  resulted  In  the  best  fit  for  extreme  value  analysis  of  low 
flows  on  most  Canadian  rivers.  However,  for  samples  with  large 
negative  skewness,  the  3-Parameter  Log  Normal  Distribution  has  pro- 
ven to  give  adequate  results  and  was  adopted  for  these  stations 
(Condie,  1983).  A  technical  discussion  of  the  fitting  procedure  Is 
given  1n  the  users  manual  and  briefly  in  Appendix  A. 
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Parameter  estimation  for  fitting  of  distributions  proceeded  in 
order  of  maximum  likelihood,  smallest  observed  drought  and 
moments.  The  procedure  utilized  for  analysis  of  each  station 
record  was  identified  on  summary  tables  together  with  the  results 
of  the  analyses  (see  Section  3.3.  and  Appendices  B-F). 

2.5  Flow  Duration  Analysis 

The  daily  discharge  data  was  extracted  for  each  station  and  sorted 
by  computer  in  order  to  derive  empirical  flow  duration  curves.  The 
flow  duration  curves  were  summarized  both  numerically  and  graphic- 
ally for  each  station  (See  Section  3.4  and  Appendices  B-F).  The 
tabular  summary  presents  the  results  of  the  flow  duration  computer 
output  at  1%  intervals  across  the  flow  duration  curve. 

2.6  Seasonal  Analysis 

In  addition  to  annual  computations  referred  to  above,  the  flow 
duration  curves  and  extreme  value  analyses  were  undertaken  on  a 
monthly  basis  to  provide  an  information  base  for  seasonal  water 
resources  investigations. 

In  order  to  reduce  the  amount  of  data  tabulation,  it  was  decided 
that  only  the  results  of  the  7O20  values  would  be  summarized  from 
the  extreme  value  output.  The  results  from  the  other  7-day 
duration  recurrence  intervals  are  available  as  part  of  the  back- 
ground files. 

Flow  duration  curves  were  also  tabulated  on  a  monthly  basis  for 
comparison  to  the  annual  curves.  A  direct  comparison  with  the 
plotted  annual  flow  duration  curve  is,  therefore,  possible  by  using 
the  numerical  tabular  summary.  Typical  results  are  discussed  in 
Section  3.4. 
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2.7    Maps 

Computer  plotting  was  utilized  as  a  cost-effective  way  of  graphic- 
ally summarizing  selected  flow  characteristics.  Base  maps  were 
digitized  for  each  region.  Data  overlays  were  then  prepared  in 
order  to  produce  two  maps  for  each  region,  the  first  map  sumariz- 
ing  selected  low  flow  characteristics  at  each  station,  and  the 
second  map  providing  a  corresponding  discharge  rate  per  unit  area. 

The  resulting  maps  are  discussed  briefly  in  Section  3.5  and  specif- 
ic maps  for  each  region  accompany  Appendices  B  to  F. 
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3.0  STUDY  RESULTS 

3.1  General 

This  section  sumnarizes  general  study  findings  and  refers  to  spec- 
ific results  given  in  separate  appendices  for  each  of  the  five 
regions.  (Appendices  B  to  F). 

3.2  Data  Base  Screening 

3.2.1  Data  Bases 

For  all  regions,  a  total  of  389  stations  were  found  to  have  10 
years  of  record  with  the  station  being  active  within  the  last  5 
years. 

Other  characteristics  of  the  data  base  were  also  tabulated  in  Sec- 
tion 2  of  Appendices  B  to  F  including,  a  brief  station  description, 
the  period  of  record  at  each  location,  whether  stream  flows  in  the 
watershed  are  considered  to  be  regulated  or  natural,  and  the  drain- 
age area  of  the  watershed. 

All  of  this  information  was  extracted  from  the  Water  Survey  of 
Canada  HYDEX  computerized  data  file  maintained  by  Environment 
Canada.  Statistical  analyses  were  undertaken  over  the  period  of 
available  record  for  both  natural  and  regulated  stations.  However, 
subsequent  to  the  data  testing  and  extreme  value  analyses  it  was 
found  that  eleven  data  series  classed  as  "regulated"  included  a 
portion  of  natural  flow  data.  These  stations  were  re-analysed  with 
results  addended  in  each  Appendix  as  appropriate.  Additional 
research  Indicated  that  three  stations  not  recorded  as  mixed  (nat- 
ural and  regulated  flows  recorded  for  the  station)  in  the  HYDEX 
file  did  consist  of  mixed  periods.  An  analysis  of  these  stations; 
02GB001.  Grand  River  below  Brantford.  02GA003  Grand  River  at  Gait, 
and,  02GA015  Speed  River  below  Guelph.  were  re-analysed  using  only 
the  regulated  period  of  data.  The  period  used  was  from  1945  (the 
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year  after  the  Shand  Dam  was  built)  to  the  present  for  the  Grand 
River  Stations  and  1975  to  present  for  the  Speed  River  data  series 
(subsequent  to  the  Guelph  Dam  construction).  The  results  for  the 
mixed  702,  ''QlO  ^"'^  ''^ZO  compared  to  the  regulated  period  only  are 
summarized  as  follows: 


Station 


Mixed  Data  Regulated  Data 

7Q2     7Qio     7Q20     70?      7Qio     7020 
m-J/s    m-^/s    m-^/s    m-^/s    m>^/s    m-^/s 


02GB001 

13.308 

7.037 

5.495 

14.980 

02GA003 

5.822 

2.175 

1.679 

9.186 

02GA015 

1.080 

0.595 

0.488 

1.399 

9.275  7.713 
5.120  4.156 
0.885    0.775 

The  results  indicate  an  increase  in  the  low  flow  values  when  using 

the  regulated  period  only,  compared  to  the  analysis  of  the  mixed 

record  period.  This  may  indicate  low  flow  augmentation,  or  trend 
in  the  data  series. 

Users  referring  to  analysis  results  for  other  regulated  stations 
should  investigate  this  aspect  in  more  detail  and  analyses  for  all 
regulated  stations  used  with  care. 

3.2.2  Screening  Results 

i)  Zero  Flows 

Prior  to  undertaking  the  extreme  value  analyses  a  total  of  23  data 
sets  containing  a  significant  number  (*)  of  zero  low  flows  for  all 
durations  were  identified  and  removed  from  the  data  base  (see  Table 
1).  For  example,  for  those  stations  where  30%  or  more  of  the 
available  extreme  flow  data  set  was  comprised  of  zero  flow,  it  was 
arbitrarily  assumed  that  use  of  this  data  base  would  result  in  a 
biased  extreme  value  analysis.  In  addition,  for  short  data  sets 
(i.e.  10  years),  reduction  by  3  years  results  in  insufficient  data 
(7  years)  with  which  to  accurately  fit  the  extreme  value 
distribution. 
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A  few  stations  which  originally  appeared  to  have  10  years  of 
records  were  found  to  have  a  significant  number  of  missing  years  in 
the  detailed  screening  process.  Those  stations  with  minimal  data 
were  removed  from  the  analysis  (refer  to  Table  1). 

ii)  Arithmetic  Overflow 

Initial  applicaton  of  the  LFA  program  encountered  an  arithmetic 
overflow  for  a  number  of  stations,  which  are  identified  in  Table 
1.  The  program  was  subsequently  modified  and  the  analysis  results 
are  included  as  addendums  to  Appendices  B  to  F  where  appropriate. 

iii)  High  Outliers 

Initially  the  analysis  of  several  stations  with  very  large  low 
flows  (greater  than  1000  m^/s)  could  not  be  completed.  These  sta- 
tions are  identified  in  Table  1.  Additional  research  was  under- 
taken and  the  LFA  program  was  subsequently  modified  to  permit  fitt- 
ing of  the  extreme  value  distribution.  These  results  are  addended 
where  appropriate  to  Appendices  B  -  F  (see  Table  1). 

iv)  Statistical  Tests 
Statistical  data  analysis  tests  were  undertaken  as  outlined  in 
Section  2.3  and  described  in  more  detail  in  Appendix  A. 3.  The 
available  test  statistics  were  recently  made  available  as  part  of 
the  LFA  (Pilon,  1988)  low  flow  analysis  program.  Previous 
application  of  these  testing  procedures  to  low  flow  series  has  been 

limited  to  date. 

In  general,  it  was  found  that  a  significant  number  of  stations 
failed  the  non-parametric  tests.  Therefore,  taken  over  the  entire 
data  base,  application  of  these  tests  has  indicated  that  the 
available  data  base  of  extreme  low  flows  may  exhibit  some  trend 
with  some  possibility  of  non-random  characteristics. 

The  data  were  further  analysed  by  subdivision  of  the  available  data 
set  according  to  length  of  record  (I.e.  ^20  years  and  vZO  years) 
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and  according  to  regulation  code.  However,  it  was  found  that 
neither  the  length  of  record,  nor  the  possible  effects  of  regula- 
tion could  account  for  the  conclusions  of  the  test  results.  One 
explanation  could  be  that  the  available  record  lengths  are  too 
short  to  permit  reasonable  application  and  interpretation  of  these 
non-parametric  test  results.  A  stronger  possibility  is  that  the 
available  low  flow  data  sets  do  exhibit  trend  and  non  random  char- 
acteristics, which  could  possibly  be  attributed  to  slow  cyclic 
change  in  groundwater  levels  or  to  climatic  trends. 

Additional  testing  was  beyond  the  scope  or  the  current  investiga- 
tions. However,  further  studies  are  recommended  since  these 
results  may  call  into  question  the  basic  assumptions  underlying 
application  of  the  extreme  value  analysis  technique  for  analysis  of 
low  flow  characteristics. 

v)  Data  Base 

Extreme  value  analyses  were  undertaken  on  an  annual  basis  for  344 
stations  and  on  a  monthly  basis  for  330  stations.  Further  analysis 
was  completed  on  the  remaining  stations  with  the  revised  program  or 
by  manual  fitting  on  the  annual  data  series. 

Fourteen  stations  were  assigned  to  more  than  one  region  since  they 
may  be  representative  of  either  region.  The  fourteen  stations 
which  appear  in  more  than  one  region  are: 


Station 

South 

Southwestern 

North 

North 

I.O. 

Central 

Eastern 

West  Central 

Eastern 

Western 

02HB011 

y 

y 

02HB012 

y 

y 

02HB013 

y 

y 

02HB016 

y 

y 

02HK002 

y 

y 

02HK003 

y 

y 

02HK0Q4 

y 

y 

02HK005 

y 

y 

02H<006 

y 

y 

04HA001 

y 

y 

04JA002 

y 

y 

04JC002 

y 

y 

04JC003 

y 

y 

04JG001 

y 

y 

TABLE  2 
Surrmary  of  Statistics  for  Data  Used  1n  the  Extreme  Value  Analysis 


Day 

Standard 

Coefficient 

Minimum 

*No.  of 

Dura- 

Mean 

Deviation 

Skew 

of 

Flow 

No.  Of 

Region 

Stations 

tion 

m^/s 

m^/s 

Variation 

m^/s 

Years 

All 

344 

1 

10.89 

3.97 

.82 

.71 

4.07 

28 

344 

3 

12.50 

4.02 

.78 

.66 

5.78 

28 

344 

7 

13.80 

4.06 

.75 

.62 

6.61 

28 

344 

15 

14.32 

4.01 

.78 

.59 

7.27 

28 

341 

30 

N/A 

3.44 

.84 

.57 

N/A 

28 

Northwestern 

67 

1 

22.05 

9.26 

.78 

.76 

7.48 

31 

67 

3 

24.35 

9.46 

.72 

.72 

8.79 

31 

67 

7 

26.73 

9.58 

.57 

.68 

9.47 

31 

67 

15 

27.97 

9.82 

.65 

.65 

9.47 

31 

67 

30 

29.35 

10.05 

.55 

.57 

11.01 

31 

Northeastern 

65 

1 

13.16 

5.60 

.79 

.70 

3.94 

31 

65 

3 

16.06 

6.02 

.78 

.64 

6.25 

31 

65 

7 

17.95 

5.62 

.62 

.58 

7.68 

31 

65 

15 

19.12 

5.67 

.63 

.53 

8.87 

31 

65 

30 

20.64 

6.05 

.59 

.48 

10.49 

31 

Southwestern 

101 

1 

1.98 

.49 

.77 

.67 

1.24 

26 

&  West 

101 

3 

2.27 

.44 

.74 

.62 

1.58 

26 

Central 

101 

7 

2.44 

.44 

.77 

.58 

1.75 

26 

101 

15 

2.61 

.49 

.85 

.56 

1.88 

26 

101 

30 

2.80 

.58 

1.03 

.57 

1.97 

26 

Central 

76 

1 

1.45 

.65 

.39 

.45 

.52 

24 

76 

3 

1.65 

.66 

.42 

.42 

.63 

24 

76 

7 

1.89 

.67 

.44 

.40 

.81 

24 

76 

15 

2.14 

.73 

.51 

.41 

1.01 

24 

76 

30 

2.53 

.91 

.61 

.40 

1.18 

24 

Southeastern 

49 

1 

27.36 

7.36 

1.62 

1.05 

12.22 

28 

49 

3 

30.86 

6.99 

1.50 

.96 

18.58 

28 

49 

7 

33.53 

7.69 

1.48 

.93 

20.97 

28 

49 

15 

33.53 

6.68 

1.45 

.89 

22.51 

28 

49 

30 

N/A 

1.51 

1.57 

.90 

N/A 

28 

•  SOME  STATIONS  APPEAR  IN  MORE  THAN  ONE  REGION. 
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The  SPSS  (Statistical  Package  Social  Science,  Norusis,  1986)  was 
used  to  produce  general  statistics  of  the  data  base  including  the 
mean,  standard  deviation  and  coefficient  of  skew  of  the  available 
low  flow  samples  for  different  durations.  These  general  statistics 
are  summarized  in  Table  2. 

With  reference  to  Table  2,  it  is  evident  that  the  average  of  the 
mean  low  flow  increases  as  the  duration  of  low  flow  increases.  The 
standard  deviation  decreases  as  a  percentage  of  the  mean  as  the 
duration  increases.  The  mean  skewness  of  the  data  decreases  with 
the  increase  in  duration. 

Flow  duration  analyses  were  subsequently  undertaken  both  on  an 
annual  and  monthly  basis.  The  flow  duration  results  are  summarized 
in  Section  3.4,  and  the  results  tabulated  on  a  regional  basis  In 
Appendices  B  to  F. 

3.3    Extreme  Value  Analysis 

3.3.1  General 

Tables  summarizing  the  results  of  the  extreme  value  analysis  for 
each  region  are  given  in  Appendices  B  to  F  for  various  recurrence 
intervals  for  each  duration.  An  example  of  the  presentation  format 
is  given  as  Table  3.  The  stations  are  identified  by  the  Water  Sur- 
vey of  Canada  station  number.  The  fitting  method  for  the  extreme 
value  distribution  is  Identified  by  a  3  letter  code,  MAX,  MOM,  SOD, 
PLN  which  stand  for;  method  of  maximum  likelihood,  method  of  mom- 
ents, method  of  smallest  observed  drought,  and  the  three-parameter 
log  normal  distribution  respectively.   (See  also  Appendix  A.l  for 
additional  Information  on  the  fitting  procedures.) 

The  next  three  parameters  appearing  in  the  summary  table  are  gener- 
al statistics  of  the  data  sample  including  the  mean  n-day  duration 
flow  (m^/s).  Standard  deviation.  Skew  (G),  and  the  coefficient  of 
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variation  of  record  (C).  The  next  two  columns  of  the  tables  give 
the  number  of  years  of  record  available  and  the  minimum  observed 
average  low  flow  for  the  particular  duration  for  the  available  data 
set.  Finally,  the  estimated  flow  for  the  various  recurrence  inter- 
vals are  listed  for  each  station. 


3.3.2  Annual 


Extreme  value  frequency  curves  were  plotted  for  each  location  and 
are  available  in  Appendices  B  to  F.  An  example  graph  is  given  in 
Tigure  1  which  summarizes  the  actual  data  and  the  fitted  curve  for 
the  1,  3,  7,  15  and  30  day  durations  for  the  Trent  River  at  Glen 
Ross,  Water  Survey  of  Canada  Station  Number  02HK004. 

A  statistical  summary  of  average  unit  area  low  flow  characteristics 
was  also  undertaken  and  is  summarized  in  Table  4.  Area  average  low 
flow  rates  (1 itres/second/kilometre^)  together  with  standard  devia- 
tion for  the  5  durations  were  tabulated  for  all  the  regions  for  the 
2  to  50  year  recurrence  intervals. 

It  is  interesting  to  note  that  in  most  cases,  the  standard  devia- 
tions are  approximately  equal  to  or  greater  than  the  corresponding 
mean  flows.   Low  flows  are  consistently  higher  than  the  Provincial 
average  in  Northwestern,  Northeastern  and  Central  Region  and  con- 
sistently lower  in  the  Southwestern  and  West  Central/Southeastern 
Regions.  However,  it  should  be  recognized  that  this  finding  may  be 
due  to  the  effects  of  the  drainage  area.  For  example,  gauged 
watersheds  in  Northern  Ontario  tend  to  be  larger  than  watersheds  in 
Southern  Ontario.  Therefore,  if  the  relationship  between  watershed 
area  and  low  flow  is  non-linear  the  corresponding  statistics  may  be 
biased. 


TABLE  4 
Unit  Area  Average  Low  Flows 


Region 


Day 
Dura- 
tion 


2 

Mean  (Std. 

Deviation) 

l/s/km2 


Recurrence  Interval  (years) 


Mean  (Std. 

Deviation) 

l/s/lan2 


10 

Mean  (Std. 
Deviation) 

l/S/kJTl^ 


20 

Mean   (Std. 

Deviation) 

l/s/lcni2 


50 

Mean   (Std. 

Deviation) 

'  /  s  /  lon^ 


All 


Northwestern 


Northeastern 


Southwestern 
&  West 
Central 


Central 


Southeastern 


1 
3 
7 

15 
30 

1 
3 
7 

15 
30 

1 
3 
7 

15 
30 

1 
3 
7 

15 
30 

1 
3 

7 
15 
30 

1 
3 
7 

15 
30 


1,59 
1.74 
1.91 
2.08 
2.31 

1.76 
1.86 
1.98 
2.08 
2.20 

1.93 
2.21 
2.48 
2.71 
3.08 

1.14 
1.24 
1.37 
1.54 
1.81 

2.18 
2.35 
2.53 
2.76 

3.14 

0.97 
1.08 
1.21 
1.32 
1.13 


1.40) 
1.54) 
1.90) 
1.79) 
1.88) 

1.17) 
1.30) 
1.53) 
1.60) 
1.66) 

1.14) 
1.37) 
1.63) 
1.73) 
1.79) 

1.17) 
1.24) 
1.34) 
1.43) 
1.59) 

1.70) 

1.83) 

1.92 

2.02) 

2.17) 

1.31) 
1.49) 
1.61) 
1.65) 
1.05) 


1.10 
1.24 
1.38 
1.52 

1.58 

1.24 
1.33 

1.43 
1.50 
1.60 

1.24 
1.49 
1.72 
1.92 
2.18 

0.76 
0.85 
0.97 
1.11 
1.30 

1.63 

1.80 
1.97 
2.16 
2.44 

0.64 
0.74 
0.84 
0.92 
0.72 


1.12) 
1.25) 
1.39) 
1.49) 
1.57) 

0.92) 

1.01) 

1.20 

1.26) 

1.32) 

0.80) 
1.01) 
1.24) 
1.36) 
1.46) 

0.96) 
1.03) 
1.12) 
1.22) 
1.36) 

1.44) 
1.56) 
1.67) 
1.78) 
1.93) 

1 
1 
1 
1 
0 


.02) 
.22) 

.33) 
.38) 
.80) 


0.89 
1.02 
1.15 
1.29 
1.43 

1.02 
1.09 
1.18 
1.24 
1.33 

0.96 
1.19 
1.40 
1.58 
1.81 

0.67 
0.70 
0.81 
0.94 
1.12 

1.38 
1.54 

1 
1 
2 


1.00) 
1.12) 
1.24 
1.34) 
1.43) 

0.80) 
0.87) 
1.03) 
1.08) 
1.14) 

0.68) 
0.86) 
1.06) 
1.19) 
1.31) 

0.88) 
0.95) 
1.03) 
1.13) 
1.26) 

1.32) 
1.42) 
1.53) 
1.65) 
1.80) 

0.89) 
1.10) 
1.21) 
1.27) 
0.70) 


0.75 
0.87 
0.99 
1.12 
1.25 

0.85 
0.92 
1.00 
1.05 
1.13 

0.77 
0.97 
1.16 
1.35 
1.55 

0.52 
0.61 
0.70 
0.84 
1.00 

1.19 
1.34 
1.50 
1.68 

1.91 


0.92 
1.02 
1.14 
1.24 

1.32 

0.73 
0.78 
0.91 
0.94 
1.00 

0.61 
0.76 
0.94 
1.07 
1.20 

0.83 
0.89 
0.96 
1.06 
1.19 

1.23 

1.30 
1.42 
1.54 
1.70 

0.80 
1.02 

1.13 
1.18 
0.64 


0.52 
0.73 
0.84 
0.97 

1.09 

0.70 
0.77 
0.82 
0.87 
0.94 

0.60 
0.78 
0.95 
1.14 
1.31 

0.44 
0.52 
0.61 
0.74 
0.90 


0.35; 
0.93) 
1.03) 
1.13) 
1.22) 

0.58) 
0.71) 
0.79) 
0.32) 
0.37) 


57) 
70) 
83) 
97) 
12) 


O.'S) 
0.83) 
0.91) 
1.01) 
1.13) 

1.14) 
1.17) 
1.29) 
1.42) 
1.59) 


.72) 
,95) 
.06) 
.11) 
.60 
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3.3.3  Seasonal 

The  7Q20  values  were  also  determined  for  each  month.  The  resulting 
analyses  are  summarized  in  Appendices  B-F.  It  is  evident  from  the 
summary  tables  that  many  of  the  lowest  of  the  7O20  generally  take 
place  in  the  summer  months. 

3.4    Flow  Duration 

3.4.1  Annual 

Flow  duration  tables  and  curves  were  produced  to  summarize  the 
percentage  of  time  the  flow  was  greater  than  or  equal  to  the  given 
value.  Flow  duration  curves  for  all  stations  are  given  in  Appen- 
dices B  to  F.  An  example  table  is  given  in  Table  5.  The  first 
column  "Per"  refers  to  the  percentage  of  the  period  of  record  that 
the  tabulated  flow  was  equalled  or  exceeded  (all  flows  are  in 
m^/s).  The  annual  flow  duration  curve  values  are  listed  in  column 
2  for  the  percentage  summarized  in  column  1.  Therefore,  the  larg- 
est daily  flow  recorded  for  this  station  up  to  1986  is  702  m^/s, 
found  in  the  1st  row  at  0  percent.  More  significantly  for  this 
study  are  the  low  flows  summarized  in  the  later  section  of  the 
table  which  have  been  exceeded  90,  95  and  100  percent  of  the  time 
for  the  period  of  record. 

The  50  percent  value  (114  m^/s)  refers  to  the  median  daily  flow  for 

the  period  of  record  of  the  flow  series  and  can  be  compared  to  the 
mean  value  145  m^/s  summarized  in  the  last  row  of  the  table.  The 
annual  flow  duration  curve  is  graphically  depicted  In  Figure  2. 
which  corresponds  to  the  numerical  values  summarized  in  columns  1 
and  2. 


TABLE  5 

EXAMPLE  SUMMARY  OF  FLOW  DURATION   ANALYSIS 

SUAMT  Tlfiii  mCN  am  CUUTION  MULYSIS  0M(OO4         TI^NT  RIVO)  AT  GLEM  ROSS 

YE>«  OF  RECWD:        23  STATION  AREA:       1200  km 2 

PO)    UHUL        MAMy    FEBRUAAT      tMCH  APKIL  MAT 


JJC        JiY       AuajsT    sEPToea  xtcbct  «vq«r  :Ecaca 


0 

7Q2.000 

425.000 

661.000 

637.000 

702.000 

515.000 

328.000 

391.000 

198.000 

340.000 

385.000 

425.000 

<39  300 

1 

518.000 

383.000 

400.000 

549.000 

864.000 

490.000 

273.000 

374.000 

1S3.000 

264.000 

378.000 

413.000 

425.000 

2 

489.000 

3BS.00O 

344.000 

530.000 

818.000 

473.000 

255.000 

289.000 

128.000 

240.000 

337.000 

398.000 

381.000 

3 

459.000 

340.000 

337.000 

515.000 

584.000 

464.000 

225. OOO 

183.000 

110.000 

225.000 

323.000 

351.000 

375.000 

4 

435.000 

319.000 

323.000 

501.000 

558.000 

459.000 

220.000 

171.000 

108.000 

198.000 

309.000 

337.000 

364.000 

5 

416.000 

294.000 

270.000 

491.000 

550.000 

439.000 

r3.ooo 

164.000 

98.100 

189.000 

272.000 

335.000 

337  CDO 

8 

396.000 

270.000 

253.000 

477.000 

534.000 

434.000 

20S.0OO 

161.000 

98.300 

167.000 

213.000 

311.000 

318.000 

7 

379.000 

266.000 

244.000 

456.000 

526.000 

431.000 

198.000 

156.000 

33.400 

150.000 

203.000 

279.000 

306  aoo 

8 

357.000 

259.000 

234.000 

442.000 

521.000 

427.000 

191.000 

152.000 

90.800 

135.000 

197. ODO 

281.000 

300  000 

9 

337.000 

253.000 

230.000 

425.000 

516.000 

416.000 

187.000 

148.000 

87.300 

129.000 

189.000 

252.000 

239. OOO 

10 

326.000 

250.000 

229.000 

419.000 

513.000 

408.000 

183.000 

138.000 

84.800 

126.000 

181.000 

249.000 

279.000 

11 

306.000 

246.000 

226.000 

410.000 

510.000 

339.000 

181.000 

127.000 

81.300 

117.000 

175.000 

246.000 

J74.QOO 

12 

290.000 

243.000 

222.000 

406.000 

504. COO 

382.000 

179.000 

121. ODO 

77.500 

105.000 

187.000 

241.000 

272.000 

13 

272.000 

240.000 

217.000 

397.000 

499.000 

376.000 

174.000 

112.000 

74.200 

33.400 

159.000 

238.000 

77TJ.0OO 

14 

263.000 

238.000 

213.000 

388.000 

496.000 

367.000 

171.000 

108.000 

70.700 

91.200 

154.000 

233. ora 

268.000 

15 

255.000 

234.000 

209.000 

374.000 

498.000 

354.000 

168.000 

103.000 

68.500 

88.800 

148.000 

229.0UJ 

266.000 

18 

247.000 

233.000 

207.000 

357.000 

490.000 

340.000 

164.000 

99.100 

68.000 

88.400 

148.000 

223. QOO 

285.000 

17 

237.000 

232.000 

203.000 

351.000 

487.000 

334.000 

161.000 

98.300 

64.000 

83.500 

140.000 

219.000 

2S2.QOO 

18 

230.000 

228.000 

202.000 

347.000 

481. QOO 

328.000 

156.000 

92.000 

61.800 

81.000 

135.000 

215.000 

260.000 

19 

224.000 

227.000 

198.000 

341.000 

475.000 

323.000 

153.000 

87.000 

60.900 

78.200 

133.000 

211.000 

258.000 

20 

219.000 

224.000 

195.000 

331.000 

470.000 

300.000 

149.000 

81.600 

59.700 

75.900 

130.000 

208.000 

254.000 

21 

213.000 

222.000 

190.000 

321.000 

487.000 

292.000 

142.000 

78.200 

58.000 

74.800 

127.000 

205.000 

252.000 

22 

208.000 

221.000 

187.000 

318.000 

467.000 

283.0OO 

138.000 

75.300 

56.900 

73.600 

124.000 

201.000 

246.000 

23 

202.000 

218.000 

184.000 

311.000 

484.000 

270.000 

133.000 

73.000 

55.400 

72.300 

121.000 

198.000 

237.000 

24 

198.000 

216.000 

183.000 

303.000 

460.000 

264.000 

130.000 

70.500 

53.800 

69.700 

119.000 

133.000 

227.000 

25 

133.000 

211.000 

181.000 

297.000 

458.000 

280.000 

128.000 

69.400 

52.800 

67.100 

115.000 

189.000 

223.000 

28 

187.000 

208.000 

180.000 

289.000 

458.000 

250.000 

120.000 

65.700 

50.700 

64.900 

113.000 

187.000 

221.000 

27 

183.000 

204.000 

178.000 

283.000 

453.000 

242.000 

116.000 

61.700 

50.000 

63.400 

111.000 

184.000 

217.000 

28 

180.000 

201.000 

175.000 

273.000 

450.000 

235.000 

112.Q00 

59.400 

47.800 

81.200 

109.000 

181.000 

215.000 

29 

177.000 

198.000 

174.000 

268.000 

447.000 

228.000 

111.000 

58.000 

47.000 

58.900 

108.000 

179.000 

213.000 

30 

173.000 

133.000 

171.000 

265.000 

445.000 

224.000 

109.000 

58.900 

46.400 

58.600 

104.000 

178.000 

212.000 

31 

169.000 

189.000 

170.000 

262.000 

439.000 

221.000 

108.000 

54.900 

45.800 

57.500 

102.000 

177.000 

206.000 

32 

165.000 

188.000 

169.000 

258.000 

438.000 

217.000 

104.000 

53.800 

44.500 

56.900 

101.000 

174.000 

204.000 

33 

162.000 

183.000 

168.000 

247.000 

430.000 

215.000 

102.000 

52.400 

44.100 

56.600 

99.800 

172.000 

202.000 

34 

159.000 

182.000 

165.000 

237.000 

429.000 

213.000 

39.100 

50.700 

43.Q00 

55.200 

98.500 

170.000 

200.000 

35 

155.000 

178.000 

1S3.000 

230.000 

422.000 

210.000 

95.700 

49.900 

42.500 

54.400 

97.400 

168.000 

199.000 

38 

152.000 

177. ODO 

162.000 

227.000 

421.000 

208.000 

91.800 

48.100 

42.000 

54.100 

96.300 

168. cm 

194.000 

37 

149.000 

174.000 

160.000 

225.000 

417.000 

207.000 

90.000 

46.400 

41.600 

53.500 

94.000 

163.000 

192.000 

38 

146.000 

172.000 

159.000 

219.000 

411.000 

204.000 

88.100 

45.400 

40.000 

53.000 

91.200 

159.000 

190.000 

39 

143.000 

170.000 

157.000 

212.000 

408.000 

201.000 

87.200 

44.500 

39.400 

51.800 

89.200 

157.000 

189.000 

40 

139.000 

168.000 

158.000 

208.000 

405.000 

137.000 

85.000 

43.800 

38.800 

51.500 

36.900 

154.000 

186.000 

41 

138.000 

168.000 

153.000 

2D3.000 

398.000 

195.000 

84.000 

43.200 

38.200 

51.000 

85.000 

152.000 

183.000 

42 

133.000 

163.000 

150.000 

198.000 

398.000 

133.000 

81.600 

42.200 

37.200 

50.100 

84.100 

150.000 

182.000 

43 

130.000 

161.000 

148.000 

133.000 

391.000 

189.000 

80.500 

41.300 

38.800 

49.800 

83.500 

147.000 

130.000 

44 

128.000 

159.000 

147.000 

189.000 

388.000 

184.000 

79.400 

40.500 

38.300 

49.200 

82.400 

146.000 

178.000 

45 

126.000 

155.000 

144.000 

188.000 

383.000 

181.000 

77.600 

39.600 

35.700 

48.700 

81.700 

143.000 

177  ODO 

46 

123.000 

153.000 

144.000 

183.000 

381.000 

178.000 

77.300 

39.100 

35.400 

47.700 

80.4CD 

142.000 

175,000 

47 

121.000 

151.000 

142.000 

179.000 

377.000 

173.000 

78.200 

38.200 

34.800 

46.400 

79.300 

137.000 

173.000 

48 

118.000 

149.000 

140.000 

189.000 

372.000 

171.000 

74.300 

37.7tn 

34.300 

45.700 

78.200 

132.000 

170,000 

49 

116.000 

148.000 

139.000 

165.000 

365.000 

170.000 

73.100 

37.100 

33.700 

44.300 

75.900 

127.000 

ISS.OOO 

TABLES  (con  t) 

SJAMY  TABLE  rnOM  an  OJUTION  AMAirSIS  OMUXM         TRENT  RIVER  AT  GLEN  WSS 

TEARS  Of  RECOC:        23  STATION  AflEA:       12000  km2 

PO)    MNJAL  JAHJAAY  FEBRUAAT      MARCH  AfWIL           MAY  JUC  JULY  tiXlST  SETTBfiCR  XTQBER  >CVQfiO)  lECaCBI 

50  114.000  146.000      138.000      162.000  357.000  187.000  71.300  38.000  33.100  44.100  74.500  '29.300  '53.000 

51  111.000  143.000      135.000      158.000  349.000  164.000  88.800  35.300  32.700  43.300  73.800  123.000  '82X0 

52  109.000  141.000      134.000      154.000  345.000  182.000  88.000  34.500  32.300  43.000  71.400  121.000  '61  000 

53  106.000  140.000      133.000      152.000  340.000  180. OX  85.700  33.700  32.300  42.500  59.300  118.000  157.3IJ 

54  104.000  139.000      132.000      145.000  337.000  158.000  64.000  33.100  31.900  41.900  67.100  116.000  '55.000 

55  103.000  138.000      131.000      142.000  331.000  156.000  82.300  32.800  31.700  41.300  88.000  115.000  153.000 

56  100. OOO  138.000      130.000      137.000  326.000  154.000  80.900  32.300  31.100  41.100  85.100  114.000  152.000 

57  97.700  133.000      129.000      138.000  315.000  151.000  59.700  32.000  -31.100  40.500  63.100  111. OOO  '51.300 

58  9S.7QO  133.000      128.000      13S.0OO  306.000  149.000  57.200  31.700  30.800  39.900  61.700  109.000  149.300 

59  33.000  131.000      127.000      132.000  303.000  148.000  55.500  31.100  30.800  38.500  60.300  107.000  '47. XO 
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3.4.2  Seasonal 

The  flow  duration  values  were  also  determined  and  tabulated  for 
monthly  flows.  Appropriate  summaries  for  all  stations  are  given  in 
Appendices  B  to  F.  This  is  also  shown  in  Table  5  in  the  example 
for  Station  02HK004,  Trent  River  at  Glen  Ross. 

The  seasonal  duration  curves  were  not  graphically  produced  for  each 
station.  However  an  example  plot  for  Station  02HK004  demonstrates 
how  the  tabulated  values  could  be  used  to  plot  monthly  curves  and 
provide  a  comparison  to  the  annual  curve  (see  Figures  3,  4  and  5). 

3.5    Maps 

Selected  low  flow  characteristics  were  summarized  on  maps  for  each 
region.  The  stations  are  located  at  the  point  of  discharge  meas- 
urement and  the  selected  data  is  summarized  in  an  information  box, 
the  format  of  which  is  depicted  in  Figure  6. 

The  stations  are  identified  by  the  7  digit  Water  Survey  Station 
Number  and  this  is  followed  by  a  regulation  code.  A  code  "R"  indi- 
cates that  data  collected  at  the  station  is  affected  by  regulation; 
the  code  "N"  means  the  station  data  are  natural  or  non-regulated. 
The  symbols  702,  ^tc.  refer  to  the  average  minimum  consecutive 
7-day  flow  (m^/s)  with  a  recurrence  interval  of  2,  5,  10  and  20 
years,  followed  by  the  minimum  one  day  flow  and  the  period  of 
record  for  the  station.  The  values  shown  on  the  right  are  the 
flows  (m^/s)  equalled  or  exceeded  for  the  available  period  of 
record  5,  50,  75,  95,  and  99  percent  of  the  time.  The  symbol  AREA 
refers  to  the  station  drainage  area  in  km^. 

Station  names  are  listed  along  with  the  station  numbers  for  identi- 
fication purposes,  on  the  map  for  each  region. 
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A  second  map  summarizes  the  above  noted  low  flow  characteristics 
expressed  as  l/s/km^.  it  may  be  possible  to  refer  to  watersheds 
with  similar  unit  runoff  rates  as  a  means  of  providing  preliminary 
low  flow  estimates  for  ungauged  watersheds.  However,  the  limita- 
tions of  area  proration  should  be  recognized.  For  example,  a  pre- 
liminary Low  Flow  Regional ization  Study  for  Southwestern  and  West 
Central  Ontario  Region  (Gumming  Cockburn  Limited,  1988)  has  found 
that  other  significant  watershed  parameters  enter  into  the  deter- 
mination of  low  flows.  Additional  investigations  are  required  in 
order  to  refine  estimating  techniques  for  ungauged  watersheds. 
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A  second  map  was  also  produced  for  each  region  which  summariz- 
ed the  above  low  flow  characteristics  expressed  as  l/s/km^. 

7.  Users  referring  to  the  analysis  results  for  regulated  stations 
should  investigate  the  effects  of  regulation  on  low  flows  in 
more  detail . 


4.2    Recommendations 


1.  Further  investigation  should  be  undertaken  to  confirm  the  app- 
licability of  available  non-parametric  tests  for  low  flow 
series. 

2.  The  possibility  of  effects  of  cyclic  changes  in  groundwater 
regime  or  climatic  changes  on  low  flows  should  be  examined  in 
future  investigations. 

3.  Additional  investigations  are  required  in  order  to  refine  low 
flow  estimating  techniques  for  ungauged  watersheds  for  each 
region. 

4.  Maps  summarizing  monthly  7O20  low  flow  characteristics  should 
be  produced  for  each  region.  This  would  be  useful  when  under- 
taking seasonal  low  flow  analyses  to  facilitate  seasonal  anal- 
yses. 

5.  Discharge  data  are  collected  continuously  by  the  Water  Survey 
of  Canada  at  each  station.  Data  analysis  and  management  tech- 
nique are  now  available  which  would  allow  efficient  updating 
of  the  present  analyses  on  a  frequent  basis.   In  our  opinion, 
the  low  flow  analyses  should  be  updated  every  three  years  in 
order  to  provide  reasonably  accurate  information  for  investi- 
gations requiring  low  flow  information. 
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4.0  CONCLUSIONS  AND  RECOMMENDATIONS 

4.1  Conclusions 

1.  The  analysis  of  non-parametric  test  results  indicate  that  the 
available  data  base  of  extreme  low  flows  may  exhibit  some 
trend  and  dependence  with  some  possibility  of  non-random  char- 
acteristics. Therefore,  some  degree  of  caution  should  be  used 
when  applying  the  results  of  this  study. 

2.  The  Gumbel  Extreme  Value  Distribution  was  found  to  adequately 
fit  the  majority  of  available  low  flow  series  for  various  low 
flow  durations.  However,  for  a  number  of  samples  with  large 
negative  skewness,  the  3  PLN  Distribution  was  adopted. 

3.  Flow  duration  analyses  were  successfully  undertaken  both  on  an 
annual  and  monthly  basis. 

4.  Extreme  value  analysis  were  undertaken  on  an  annual  basis  for 
344  stations  and  on  a  monthly  basis  for  330  stations  using  the 
LFA  program.  Graphical  analyses  were  undertaken  for  the 
remainder  of  the  stations. 

5.  It  was  found  that  area  average  low  flows  are  consistently 
higher  than  the  Provincial  average  in  the  Northwestern, 
Northeastern  and  Central  Regions  and  consistently  lower  in  the 
Southwestern/West  Central  and  Southeastern  Regions. 

6.  A  map  summarizing  the  following  low  flow  characteristics  vas 
produced  for  each  region: 

-  7  day  average  extreme  values  for  the  2,  5,  10  and  20  year 
recurrence  intervals 

-  flow  duration  values  which  were  equalled  or  exceeded  over 
the  available  period  5,  50,  75,  95  and  99  percent  of  the 
time 
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annual  and  monthly  basis. 

4.  Extreme  value  analysis  were  undertaken  on  an  annual  basis  for 
344  stations  and  on  a  monthly  basis  for  330  stations  using  the 
LFA  program.  Graphical  analyses  were  undertaken  for  the 
remainder  of  the  stations. 

5.  It  was  found  that  area  average  low  flows  are  consistently 
higher  than  the  Provincial  average  in  the  Northwestern, 
Northeastern  and  Central  Regions  and  consistently  lower  in  the 
Southwestern/West  Central  and  Southeastern  Regions. 

6.  A  map  summarizing  the  following  low  flow  characteristics  -as 
produced  for  each  region: 

-  7  day  average  extreme  values  for  the  2,  5,  10  and  20  year 
recurrence  Intervals 

-  flow  duration  values  which  were  equalled  or  exceeded  over 
the  available  period  5,  50.  75,  95  and  99  percent  of  the 
time 
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A  second  map  was  also  produced  for  each  region  which  surmariz- 
ed  the  above  low  flow  characteristics  expressed  as  l/s/km^. 

7.  Users  referring  to  the  analysis  results  for  regulated  stations 
should  investigate  the  effects  of  regulation  on  low  flows  in 
more  detail. 


4.2    Recommendations 


1.  Further  investigation  should  be  undertaken  to  confirm  the  app- 
licability of  available  non-parametric  tests  for  low  flow 
series. 

2.  The  possibility  of  effects  of  cyclic  changes  in  groundwater 
regime  or  climatic  changes  on  low  flows  should  be  examined  in 
future  investigations. 

3.  Additional  investigations  are  required  in  order  to  refine  low 
flow  estimating  techniques  for  ungauged  watersheds  for  each 
region. 

4.  Maps  summarizing  monthly  7020  ^o^  ^^°^  characteristics  should 
be  produced  for  each  region.  This  would  be  useful  when  under- 
taking seasonal  low  flow  analyses  to  facilitate  seasonal  anal- 
yses. 

5.  Discharge  data  are  collected  continuously  by  the  Water  Survey 
of  Canada  at  each  station.  Data  analysis  and  management  tech- 
nique are  now  available  which  would  allow  efficient  updating 
of  the  present  analyses  on  a  frequent  basis.   In  our  opinion, 
the  low  flow  analyses  should  be  updated  every  three  years  in 
order  to  provide  reasonably  accurate  information  for  investi- 
gations requiring  low  flow  information. 
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Introduction 

Statistical  distributions  allow  estimates  of  probability  of  exceedance  of 
events  to  be  made  by  analytical  techniques.  Various  methods  are  avail- 
able, and  a  good  discussion  is  given  by  McMahon  and  Mein  (ref  9). 

At  a  gauged  location,  average  low  flows  can  be  determined  for  selected 
durations  of  time  for  each  year  of  record.  For  example,  the  corresponding 
minimum  average  consecutive  15-day  low  flow  can  be  determined  for  each 
annual  record.  The  series  of  15-day  average  annual  low  flows  can  then  be 
ranked  lowest  to  highest  and  an  extreme  value  analysis  undertaken  using  a 
theoretical  distribution. 

The  Gumbel  Extreme  Value  Distribution  is  commonly  used  for  fitting  low 
flow  frequency  curves  to  the  available  data. 

The  Gumbel  Type  III  Distribution 

This  is  a  variation  of  Fisher  and  Tippett's  third  asymnptotic  distribution 
of  extreme  values  and  is  sometimes  referred  to  as  the  Weibull  distribu- 
tion. In  view  of  his  many  contributions  it  is  often  referrred  to  in 
hydrology  as  the  Gumbell  II.  The  probability  density  function  is: 


4>(x)  =   a 

u  -  e 


X  -  e 

u  -  6 


exp 


X  -  e 

u  -  e 


(1) 


where  e  is  the  lower  bound  parameter 
u  is  the  characteristic  drought 
and   a  is  the  shape  parameter 


Various  methods  are  available  for  determining  the  distribution  parameters 
for  a  particular  data  set.  See  references  2  and  6  for  example. 


The  Probability  Function 


The  density  function  is  integrable  and  gives  the  distribution  function 


F  (x)  =  1  -  exp 


X  -  e 

U  -  6 


(2) 


which  gives  the  probability  of  non-exceedance  of  x. 

Since  it  is  more  common  to  require  x  for  a  given  probability  of 

non-exceedance,  a  simple  re-arrangement  gives 


X  =  e  +  (u  -  e)  (  -  1 


n  fl  -  F(x)  ~| 


1/a 


(3) 


There  are  occasions  when  the  sample  series  of  low  flows  can  have  a  large 
negative  coefficient  of  skewhess,  and  since  the  Gumbel  Type  III  distribu- 
tion cannot  have  a  skewness  of  less  than  -1.14,  then  the  distribution 
cannot  be  fitted. 

There  are  insuffucient  natural  samples  available  with  such  low  skewness  to 
make  a  firm  choice  of  an  alternative  treatment,  but  from  the  few  avail- 
able, it  was  found  that  the  negatively  skewed  three-parameter  lognormal 
distribution  provided  an  acceptable  alternative. 


If  y  =  In  (a-x)  is  normally  distributed,  then  the  probability  density 
function  of  x  is  given  by 

2  ^ 


(t.(x)  = 


1 


exp 


(a-x)o  yzn 


j  -_L  [in  (a-x)  -tnj 


(4) 


where  m  and  o  are  respectively  the  mean  and  standard  deviation  of  the 
series  In  (a-x)  and  a  Is  an  lower  boundary  parameter.  This  form  of  the 
distribution  can  only  have  negative  skewness.  The  equation  1n  the  form 
In  (x-a)  Is  for  distributions  with  positive  skewness. 


Taking  moments,  re-arranging  and  replacing  them  by  their  sample  estimates 
gives: 

k^  +  3k  -  gi  =  0  (36a)   (5) 

and  after  solving  for  k 

a  =  X  -  s/k  (37a)   (6) 

m  =  ln(-s/k)  -  l/21n(k2+l)  (38a)   (7) 

o2  =  ln(k2+l)  (39a)   (8) 

Thus,  the  distribution  is  completely  defined  and  the  T-year  low  flows  can 
be  computed  from: 

Qt  =  a  -  exp(m  ot)  (9) 

Additional  details  are  discussed  in  the  user's  manual. 


Source:  Condie,  R.,  L.  C.  Cheng 

"Low  Flow  Frequency  Analysis  -  Program  LOFLOW", 
Water  Resources  Branch,  Inland  Waters  Directorate 
Environment  Canada,  Ottawa,  1987 
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APPENDIX  A. 2 
FLOW  DURATION  ANALYSIS 


A  flow  duration  curve  is  a  plot  of  the  flow  of  a  stream  against  the  per- 
cent of  time  the  indicated  flow  was  equalled  or  exceeded  during  the  period 
covered  by  the  available  flow  data.  This  curve  is  extremely  useful  for 
hydro  power  studies  and  for  characterizing  the  local  streamflow  regime. 

An  empirical  procedure  was  used  to  analyse  each  sample  discharge  record. 
Flow  duration  curves  are  derived  by  rearranging  the  available  daily  or 
monthly  flow  data  in  order  of  magnitude.  The  total  time  period  represents 
1007.  of  the  time.  Therefore,  by  definition  of  the  procedure  used  here, 
the  largest  value  is  exceeded  0  percent  of  the  time  and  the  smallest  value 
is  exceeded  1007.  of  the  time. 

The  flow  duration  curve  represents  the  percent  of  time  that  a  specific 
discharge  occurs  at  that  location.  However,  the  curve  does  not  indicate 
the  period  of  time  in  the  year  when  the  flow  is  less  than  or  equal  to  the 
selected  value.  Therefore,  in  some  instances,  it  is  also  useful  to  devel- 
op flow  duration  curves  on  a  monthly  basis,  that  is,  all  data  for  the 
month  of  January  over  the  entire  period  of  record  is  analysed  independent- 
ly to  produce  a  flow  duration  curve  representative  of  flow  conditions  in 
January. 

Single  station  flow  duration  analyses  and  corresponding  monthly  flow  dura- 
tion curves  have  been  determined  at  all  relevant  stream  and  flow  locations 
across  the  Province. 
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APPENDIX  A. 3 

NONPARAMETRIC  TESTS  FOR  INDEPENDENCE, 

TREND  AND  RANDOMNESS 


A. 3.1    Test  Description 


This  appendix  briefly  summarizes  the  functions  evaluated  in  the 
package  and  gives  the  methods  used  to  determine  their  statistical 
significance.  Statistical  tables  are  provided  for  ease  of 
reference. 


A. 3. 1.1  Introduction 


Any  statistical  test  of  significance  will  generally  be  made  using 
the  following  steps: 

a)  State  the  null  hyposthesis,  Hq.  For  instance  in  split  sample 
tests,  the  null  hypothesis  may  be  that  there  is  no  difference 
between  the  sample  means. 

b)  Choose  a  significance  level,  a. 

c)  Choose  an  appropriate  statistical  test.  In  this  program  all 
tests  are  nonparametric. 

d)  Compute  the  test  statistic. 

e)  The  sampling  distribution  of  the  test  statistic  is  known  and 
has  been  tabulated,  and  the  chosen  significance  level  then 
defines  the  region  of  rejection. 

f)  If  the  computed  test  statistic  lies  in  the  region  of 
rejection,  then  the  null  hypothesis  Is  rejected. 

Consider  now  the  four  tests  in  this  program. 


A. 3. 1.2  The  Spearman  Rank  Order  Serial  Correlation  Coefficient  for 
Independence 

If  the  series  Q.  with  i  ranging  from  1  to  N  is  put  in 
chronological  order,  two  time  series  are  formed  and  their 
respective  ranks  computed: 

Oi,  02  On-1  by  X.,  the  rank  of  0^;  i  =  l;  N-1 

and  O2.  Q3  On  by  y. ,  the  rank  of  0.;  i=2,  N 

then  Spearman  rank  order  serial  correlation  coefficient  is 

S^=  I   (Zx2+  Zy2  .  2d2)  (2x2  2^2)4  (Ij 

where  Zx?  =  'T\'^-m)l\2-   2Tx 

Sy^  =  (m3-m)/12-  ZT 

d.  is  the  difference  in  rank  between  x.  and  y. 
m  =  N-1 


and  the  summations  are  over  the  m  pairs  of  x.,  y.. 


Ignoring  for  the  moment  the  terms  in  T  and  putting  them  at  zero, 
equation  (a)  becomes 

S^=  l-(6Zd?)  (m^-m)  (2) 


the  more  familiar  form  of  the  Spearman  rank  correlation 
coefficient. 

The  terms  in  T  adjust  for  tied  ranks  and  are  computed  as 
follows.  If  for  instance  three  observations  in  the  x  series  were 
tied  for  ranks  17,  18,  and  19  then  each  observation  is  given  the 
rank  18;  if  two  were  tied  for  ranks  24  and  25,  then  each  is 
ranked  24.5. 


For  each  tied  set,  T  is  computed  from 


T^  =  (J3  -  J)/12 


where  J  is  the  number  of  observations  tied  at  a  give  rank. 

ZT  and  2T  are  defined  by  the  extension  of  the  foregoing. 
X      y 

When  N  is  10  or  greater,  then  the  function 

t  =  S^  [{m-2)/(l-s2)]'^  (3) 

is  distributed  like  Student's  t  with  m-2  degrees  of  freedom.  A 
one-tail  test  must  be  used. 

A. 3. 1.3  The  Spearman  Rank  Order  Correlation  Coefficient  Test  for  Trend 

If  the  series  0^  with  i  ranging  from  1  to  N  is  put  in 
chronological  order,  ranks  assigned  and  denoting  the  series 

Ql,  02  On  by  y^,  the  rank  of  0^ 

and  1,  2,  N  by  x. ,  the  sequential  order  of  0^ 

then  the  Spearman  rank  order  correlation  coefficient  r  is 
calculated  as  in  equation  1,  except  that  m  =  N,  T  =  0,  and  the 

summations  are  taken  over  the  N  pairs  of  x. ,  y,. 

For  N  =  10  or  greater,  then  the  function 

t  =  r^  [(N-2)/(l-r2)]^  (4) 

is  distributed  like  Student's  t  with  N-2  degrees  of  freedom.  The 
null  hypothesis  Is  that  there  Is  no  trend,  either  upward  or 
downward  with  time,  and  so  a  two-tall  test  Is  used. 


A. 3. 1.4  Runs  Above  and  Below  the  Median  for  General  Randomness 

This  randomness  test  is  based  on  the  order  or  sequence  in  which 
the  Individual  scores  or  observations  were  obtained.  Actually, 
the  test  is  based  on  the  number  of  runs  which  a  sample  exhibits. 

A  run  is  defined  as  a  succession  of  identical  symbols  which  are 
followed  and  preceded  by  different  symbols  or  by  no  symbols  at 
all. 

The  total  number  of  runs  in  a  sample  of  any  given  size  gives  an 
indication  of  whether  or  not  the  sample  is  random.  If  very  few 
runs  occur,  a  time  trend  or  some  bunching  due  to  lack  of  indepen- 
dence is  suggested.  If  a  great  many  runs  occur,  systematic 
short-period  cyclical  fluctuations  seem  to  be  influencing  the 
sample. 

For  example,  once  the  median  of  the  sample  has  been  determined, 
each  observation  can  be  labelled  as  being  above  and  equal  to  or 
below  and  equal  to  the  median.  If  "A"  represents  above  and  equal 
to  the  median  and  "B"  represents  below  and  equal  to  the  median, 
then  a  sample  may  be  ordered  as 

AABBBABBBBAABA 

(A  run  represents  a  succession  of  identical  symbols).  For  our 
example,  each  run  is  underscored  and  numbered  consecutively: 


AA 

BBB 

A 

BBBB 

AA 

B 

A 

1 

2 

3 

4 

5 

6 

7 

This  sample  begins  with  2  observations  above  or  equal  to  the 
median,  followed  by  a  run  of  3  observations  below  or  equal  to  the 
median,  etc. 

Seven  runs  are  observed  in  all:  that  is,  the  total  number  of 
runs  above  and  below  the  median  RUNAB,  is  7.  If  n]  represents 


the  number  of  events  of  type  A,  then  n^  =  6.  If  n2  denotes  the 
number  below  the  median,  type  B,  then  n2  =  8.  Thus,  the  number 
of  observations  is  equal  to  (n^  +  n2). 

In  order  to  apply  this  run  test,  one  must  determine  nj,  n2  and 
RUNAB. 


The  null  hypothesis,  H  ,  is  that  the  A's  and  B's  occur  in  random 
order.  The  alternate  hypothesis,  H.,  is  that  the  order  of  the 
A's  and  B's  deviates  from  randomness. 


When  either  n^  or  n2  is  greater  than  20,  the  sampling  distribu- 
tion of  RUNAB  tends  to  normality  with 


IrUNAB  -  [(2nin2)/(ni+n2)+l]  | 

(5) 


2        ^i 
{2nin2(2nin2-ni-n2)/[(ni+n2)  (ni+n2-l)]^ 


where  z  is  an  N(0,1)  variate  as  described  in  Table  A. 4.  This 
package  uses  a  region  of  rejection  defined  by 

z  greater  than  1.96  for  a  =  .05 

z  greater  than  2.575  for  a  =  .01 


A. 3. 2    Example  of  Tests 


The  non-parametric  test  referred  to  in  Section  2.3  were  performed 
on  all  extreme  value  data  sets.  The  data  for  the  7-day  duration 
for  Station  02BD002  which  "failed"  the  non-parametric  tests  is 
tabulated  in  Table  1  and  comparatively  the  data  for  Station 
02AB009  which  passed  all  the  non-parametric  tests  is  given  in 
Table  2.  The  results  of  the  tests  for  independence,  trend  and 
randomness  for  both  stations  are  given  in  Tables  3  and  4  respect- 
ively. 


Source:  P1lon,  P.J.,  R.  Condle.  K.  D.  Harvey 

"Consolidated  Frequency  Analysis  Program  -  CFA".  July  19S5 
Water  Resources  Branch,  Inland  Waters  Directorate. 
Environment  Canada  -  Ottawa.  1985 
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TABLE  2 

7   DAY   DURATION   FLOWS  USED 
FOR  EXTREME   VALUE  ANALYSIS 
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0260002 


7  DAY  LDU  FLQU  I1EAH  OISCH  IH  PERIOD  JAH  01  DEC  31 


vflftTIWC 

7  DAY 

AiCEHDIHG 

CUnULAT 

RETiJRH 

tlQWTH 

YEAR 

flEAH  FLDU 

ORDER 

RAMK 

PROBABIL. 

PERIOD 

fcns) 

(cns) 

•  i-ESRS; 

10 

1924 

17  0000 

5  3860 

1 

1.01 

98  47 

9 

1925 

20  9710 

9  6090 

2 

2.70 

37  00 

A 

1926 

19  0140 

15  2710 

3 

4.39 

22  77 

9 

1927 

23  6430 

15.9710 

4 

6.08 

16  44 

3 

1928 

2^  7000 

17  0000 

5 

7  77 

l2  87 

9 

1933 

15  2710 

19  0140 

4 

9  46 

10  57 

1 

1934 

21.3140 

20  6430 

7 

11  15 

6  97 

9 

1935 

27   3710 

20  9710 

3 

12  84 

7  79 

12 

1936 

2'j  3140 

21  1570 

9 

14  53 

6  88 

4 

1937 

2  3  1000 

21  2140 

10 

U.22 

4  17 

10 

1938 

23  0570 

21  3140 

11 

17  91 

5  58 

8 

1939 

25  1860 

22.7710 

12 

19.59 

5  10 

3 

1940 

25  3000 

23.0570 

13 

21.28 

4  70 

8 

1941 

24  2430 

23.1000 

14 

22.97 

4  35 

7 

1942 

27  5570 

23.5000 

15 

24.66 

4  05 

12 

1943 

25  9430 

23  8860 

16 

26  35 

3  79 

1 

1944 

23.8860 

24  2430 

17 

28  04 

3  57 

3 

1945 

24.5430 

24.5430 

18 

29  73 

3  36 

10 

1946 

26  1570 

25  1860 

19 

31.42 

3  18 

12 

1947 

21  1570 

25.3140 

20 

33.11 

3  02 

10 

1948 

20  6430 

25  6430 

21 

34.80 

2  87 

12 

1949 

23  5000 

25.8000 

22 

36.49 

2  74 

2 

1950 

31.1710 

25.9430 

23 

38.18 

2.62 

9 

1951 

50.0860 

26  1570 

24 

39  86 

2  51 

11 

1952 

44.0570 

26  7000 

25 

41  55 

2  41 

10 

1953 

31  1000 

27  2000 

26 

43.24 

2  31 

11 

1954 

36.2140 

27  3710 

27 

44  93 

2  23 

9 

1955 

21.2140 

27  5570 

28 

46  62 

2  14 

4 

1956 

32.2290 

27  7000 

29 

48  31 

2  07 

7 

1957 

39  3710 

29  7000 

30 

50  00 

2  00 

8 

1958 

31.8000 

30  4430 

31 

51  69 

1  93 

6 

1959 

42  4000 

30  8710 

32 

53  38 

1  87 

10 

1960 

30  4430 

31  1000 

33 

55  07 

1  82 

7 

1961 

41  0000 

31  1710 

34 

56.76 

1  76 

8 

1962 

42  1860 

31  8000 

35 

58  45 

1  71 

11 

1963 

38  2290 

32  2290 

36 

60  14 

1  64 

6 

1964 

30.8710 

32  4570 

37 

61  82 

1  42 

9 

1965 

61  0140 

34  9000 

38 

63  51 

1  57 

10 

1966 

36  6290 

36  2140 

39 

65  20 

1  53 

12 

1967 

48  7430 

36  4290 

40 

66  89 

1  49 

3 

1968 

32  4570 

38  2290 

41 

68  58 

1  46 

9 

1969 

5  3860 

39  3710 

42 

70  27 

I  42 

12 

1970 

51  04  30 

39  6430 

43 

71  96 

1  39 

11 

1971 

34  9000 

39  7570 

44 

73  45 

1  36 

5 

1972 

45  1860 

41  0000 

45 

75  34 

I  33 

12 

1973 

58  3000 

•n  4710 

44 

77  03 

1  30 

5 

1974 

39  7570 

42  1860 

47 

78  72 

1  27 

11 

1973 

15  9710 

42  4000 

41 

80  41 

I  24 

10 

1976 

22  7710 

44  0570 

49 

82  Of 

1  22 

1 

1977 

9  4090 

45  1840 

50 

83  78 

1  19 

10 

1978 

41  4710 

44  4710 

51 

85  47 

I  17 

9 

1979 

53  1430 

48  74  30 

52 

87  16 

1  15 

6 

1980 

39  6430 

50  0840 

53 

88  85 

1  13 

3 

1981 

27  7000 

50  2290 

34 

90  54 

I  10 

"J 

1982 

26  7000 

51  0430 

53 

92  23 

I  08 

7 

1983 

55  14  30 

53  1430 

36 

93  92 

1  06 

10 

1984 

50  2290 

55  1450 

37 

93  61 

1  OS 

'j 

1985 

44  4710 

58  3000 

38 

97  30 

t  03 

8 

1986 

27  :ooo 

41  OHO 

59 

98  99 

1  01 

TABLE  3 

RESULTS  OF  STATISTICAL  TESTS 
FOR  STATION  02AB009 


---  SPEARnflN  TEST  FDR  IHOEPEHOEXCE  — 

0.'AtOu913007      SHEkiAHOQUflH  RIUER  AT  SUNSHIHE 

ftXHUAL  mmnun  OAILY  FLDU  series    38  TQ   86  ORAIHAGE  AREA  =  2800  000 

SPfHWIIflU  RAHK  ORDER  SERIAL  CORRELATION  CQEFF  =  0  068  D  F  =   26 

CORRESPONDS  TO  STUOENTi  T  =  0  347 

CRITICAL  T  UALUE  AT  5/  LE'JEL  ^  1  ?G6  NOT  SIGNIFICANT 

-   -  -     U       -       =  2  479  NOT  SIGNIFICANT 

Inrerprirt ation  The  null  hypothesis  is  that  the  sen jl(lig-one)  correlatio« 
15  zero 

At  rhe  V  level  of  significance,  the  correlation  is  not  Signif icinllg 
different  fron  zero   That  is,  the  data  do  not  display  Sigiificant  serial 


—  SPEARflAH  TEST  FDR  TREND  — 

02aE:00913007      SHEtAHDOUAM  RIUER  AT  SUMSHIME 

ANNUAL  niHinun  DAILY  FLOU  SERIES    38  10   86  DRAIHACE  AREA  =  2800  000 

SPEARHAN  RANK  ORDER  CORRELATION  COEFF  =  -0.227  OF  =   27 

CORRESPONDS  TO  STUDENTS  T  =  -1  209 

CRITICAL  T  UALUE  AT  'i/.   LEUEL  =  -2.052  NOT  SICHIFICANT 

-   -   -  1/   -  =  -2.771  NOT  SIGNIFICANT 

InterpretatioA:  The  null  h^^othesis  is  that  tke  serial(laj-ote)  correlation 
15  zero. 

At  the  5/  level  of  significance,  the  correlatioi  is  not  signif icantly 
iiffereot  fron  zero.  That  is,  the  data  do  not  display  sigatficaat  trend. 


--  RUN  TEST  FOR  GENERAL  RANDQtlNESS  — 

02A{;00913007      SHEtANODUAN  RIUER  AT  SUNSHINE 

ANNUAL  niNIflUn  DAILY  FLOU  SERIES    38  TO   86  DRAINAGE  AREA  =  2800.000 

THL  NUt1E;tR  OF  RUNS  AfcOUE  AND  HELQU  THE  HEDIAN(RUNAEi)  =  11 

THE  NUnBER  Of  OBSERUATIDNS  ABOUE  THE  f1£DIAN(Nl)  =  14 

THE  NUflPER  OF  OBSERUATIDNS  BELDU  THE  I1EDIAM(N2)  =  14 

Range  at  3/  level  of  significance  10  to  20     NOT  SIGNIFICANT 

Interpretation:  The  null  hypothesis  is  that  the  data  are  randon. 

At  the  5/  level  of  siqnif icance,  the  null  hypothesis  caniot  be 
rejected   That  is,  the  sanple  is  siynif icaatly  raidon 


TABLE  4 

RESULTS  OF  STATISTICAL  TESTS 
FOR  STATION  02AB009 


—  SPEflRllflH  TEST  FUR  IHDEPENOEHCE  — 

02t:D0021Ju07      niCHlPICQTEH  RIUER  AT  HIGH  FftLLS 

AHHUAL  niHIflUn  DftlLV  FLQU  SERIES    24  TQ   86  ORAIHAfiE  AREA  =  3130  000 

SPEftRdflN  RAKK  ORDER  SERIAL  CDRRELATIDH  COEFF  =  0  «6  0F.=   36 

CORRESPOHOS  TO  STUOEHTS  T  =  3.839 

CRITICAL  T  UALUE  AT  V.   LEUEL  =  1.674  SlgHIFICAHT 

-   -   -  1/   -   =  2  397  SICHIFICAHT 

Interpretition:  The  null  hijpothesis  is  that  tbe  seriil(la9-o»e)  correlatioi 
15  zero. 

At  the  1/  level  of  significance,  the  correlation  is  sijnif icantl? 
different  fron  zero.   That  is.  the  4ata  display  highlji  significant 
serial  dependence. 


—  SPEARttAK  TEST  FDR  TREHD  — 

02l:D00213007      I1ICHIP1CQTEH  RIUER  AT  HICH  FALLS 

ANNUAL  (liHinun  DAILY  FLQU  SERIES    24  TO   86  DRAINAGE  AREA  =  5130.000 

SPEARHAX  RANK  ORDER  CQRRELATIDH  CQEFF  =  -0  528  D.F.=   57 
CORRESPONDS  TO  STUDENTS  T  =  -4  693 

CRITICAL  T  UALUE  AT  5/  LEUEL  =  -2.003  SIGNIFICANT 

-   -   -  IX   -   =  -2  667  SIGNIFICANT 

Interpretation:  The  null  hypothesis  is  that  the  seriaKlaj-one)  correlation 
15  zero. 

At  the  1/  level  of  significance,  the  correlatioi  is  sigiif iciitlj 
different  fron  zero.   That  is  the  data  display  hijklt)  significant  trend 


---  RUN  TEST  FDR  GENERAL  RANDOttNESS  — 

021000213007      I11CH1PIC0TEN  fllUER  AT  HIGH  FALLS   ^^  ^^^^    c,  ,„  ,in« 
ANNUAL  niNinun  daily  FLQU  SERIES    24  TO   86  DRAINAGE  AREA  =  5130  000 

THE  NUnttR  OF  RUNS  ABOUE  AND  KELDU  THE  nEDlAN<RUNAb)  =  11 
THE  NUHBER  of  OBSERUATIQNS  ABQUE  the  nEDIAN(Nl)  =  29 
THE  NUnPER  OF  QBSERUATIONS  BELOU  THE  nEDIAN<N2)  =  29 

(NOTE  2  IS  THE  STANDARD  NQRHAL  UARIATE  > 

For  this  test,  Z  =    5  034 

Critical  Z  value  at  the  5/  level  =  1  960  cicillr lr2!Il 

Critical  Z  value  at  the  1/  level  --   2  575  SIGNIFICANT 

Interpretation  The  null  hypothesis  is  that  the  data  are  rindon 

At  the  1/  level  of  significance,  the  null  hgpotkeiis  cm  he  rejected 
That  IS,  the  sanple  is  not  si gnif icantlg  raidon 


A. 3. 3    Graphical  Representation 

Graphical  procedures  are  available  within  the  LFA  (Pilon, 
Jackson,  1987)  program  which  visually  depict  some  of  the  results 
of  the  non-parametric  tests  as  shown  in  figures  1  to  6. 

Figures  1  and  2  show  that  there  is  an  apparent  trend  for  Station 
02BD002  when  compared  with  the  rank  versus  time  graph  of  Figure 
2.  These  results  correspond  to  the  statistical  test  results. 
The  increase  in  low  flows  as  a  function  of  time  is  very  apparent 
for  the  data  for  Station  02BD002  (see  Figure  3)  when  compared  to 
the  data  from  Station  02AB009  shown  in  Figure  4.  Figures  5  and  6 
graphically  depict  the  probability  density  function  by  displaying 
a  histogram  of  the  number  of  occurrences  of  between  ranges  of  low 
flows.  It  can  be  seen  that  both  stations  are  positively  skewed 
for  the  7-day  duration  data. 

A. 3. 4    Summary  of  Tests  Results 

Test  results  are  summarized  by  region  for  significant  levels  of 
1%  and  5%  in  appendices  B  to  F,  Section  2.0.   It  can  be  seen 
from  the  summary  tables  (Tables  2  to  6)  that  the  data  for 
stations  in  the  Central  Region  are  random  at  the  17.  level  (refer 
to  Table  2). 

All  other  test  results  indicate  that  a  large  number  of  the  data 
sets  "fall"  the  non-parametric  test.  This  could  indicate  some 
dependence,  trend  and  non-randomness.  Hence,  some  degree  of 
caution  should  be  used  when  applying  the  results  of  this  study. 

To  examine  effects  of  length  of  record  and  regulation,  the 
analysis  was  completed  using  criteria  of  >_20  years  of  record,  <20 
years  of  record  and  regulated  vs.  non-regulated  station  data. 
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TABLE  2 
Sunrnarv  of  Data  Screening  All  Stations 


Northeastern  Region 


Inaepenaence 

T 

rend 

Sana 

mness 

Day 

11 

51 

11 

5X 

11 

•         j 

Duration 

Sig. 

Not 

Per. 

sig. 

Not 

Per. 

sig. 

Not 

Per. 

Sig. 

Not 

Per. 

sig. 

Not 

Per. 

Sig.  '    tc 

:   =•-.  1 

1 

53 

7 

0 

47 

13 

0 

5D 

10 

0 

17 

0 

46 

19 

0 

a 

3 

51 

9 

0 

44 

13 

0 

48 

12 

0 

17 

0 

47 

18 

0 

46 

; 

7 

50 

10 

0 

44 

16 

0 

50 

10 

0 

19 

0 

47 

18 

0 

48  1  ; 

0   1 

15 

51 

9 

0 

46 

16 

0 

53 

7 

0 

13 

0 

46 

19 

0 

48  1  : 

'  '  3 

30 

54 

S 

0 

46 

14 

0 

50 

10 

0 

16 

0 

49 

16 

0 

48  1  : 

D   ! 

•  TOTAL 

259 

41 

255 

75 

251 

49 

213 

87 

235 

90 

226     : 

9 

Northwestern  Region 


Independence 

T 

■end 

Sanaanness 

Day 

n 

51 

IX 

51 

11 

=• 

Duration 

S19. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Stg. 

Not 

Per. 

Sig. 

Not 

=  »-. 

1 

47 

20 

0 

42 

25 

0 

59 

8 

0 

49 

18 

0 

31 

36 

0 

46 

21 

3 

3 

48 

19 

0 

44 

23 

0 

60 

7 

0 

53 

14 

0 

29 

38 

0 

47 

20 

0 

7 

49 

18 

0 

45 

22 

D 

60 

7 

0 

53 

14 

0 

30 

37 

0 

SO 

17 

3 

15 

51 

16 

0 

43 

24 

0 

62 

5 

0 

54 

13 

0 

31 

36 

0 

51 

30 

50 

15 

0 

42 

25 

0 

62 

5 

0 

50 

17 

0 

30 

37 

0 

48 

:5 

3 

•TOTAL 

247 

88 

216 

119 

303 

32 

259 

76 

151 

184 

242     53 

Southwest  and  West  Central  Region 


Independence 

Trend 

aanoomness 

Day 

11 

51 

11 

51 

11 

51 

Duration 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

1 

72 

26 

0 

66 

32 

0 

57 

40 

0 

54 

0 

17 

84 

0 

67 

3 

3 

71 

27 

0 

68 

30 

0 

54 

44 

0 

55 

0 

18 

83 

0 

68 

0 

7 

72 

26 

0 

66 

32 

0 

58 

40 

0 

56 

0 

19 

82 

0 

71 

0 

15 

74 

24 

0 

67 

31 

0 

57 

41 

0 

54 

0 

17 

84 

0 

73 

0 

30 

71 

27 

0 

65 

33 

0 

54 

44 

0 

57 

0 

17 

84 

0 

70 

0 

.  TOTAL 

360 

130 

332 

158 

280 

210 

214 

276 

88 

417 

349 

156 

Central  Region 


Independence 

Trend 

Sinaomness 

Day 

11 

51 

11 

51 

11 

il 

Duration 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Stg. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

1 

64 

11 

0 

56 

19 

0 

53 

22 

0 

37 

0 

1 

:'5 

52 

62 

14 

0 

3 

65 

10 

0 

56 

19 

0 

52 

23 

0 

41 

0 

1 

75 

52 

60 

16 

0 

7 

64 

11 

0 

56 

19 

0 

54 

21 

0 

37 

0 

1 

75 

52 

62 

14 

0 

15 

66 

9 

0 

58 

17 

0 

49 

26 

0 

39 

0 

0 

76 

100 

60 

16 

0 

30 

67 

a 

0 

56 

19 

0 

49 

26 

0 

38 

0 

1 

75 

52 

66 

10 

0 

•  TOTAL 

326 

49 

282 

93 

257 

118 

192 

183 

4 

376 

310 

'0 

Southeastern  Region 


Independence 

Trend 

■landoantst 

Day 

11 

51 

11 

51 

11 

SI 

Duration 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

SI9. 

Not 

Per. 

SI9. 

Not 

Per. 

S'9. 

Not 

Hr. 

1 

39 

0 

9 

0 

0 

36 

10 

0 

22 

27 

0 

N 

» 

3 

39 

0 

10 

0 

0 

3S 

a 

0 

21 

27 

0 

V 

\t 

7 

40 

0 

7 

0 

0 

41 

s 

0 

23 

2( 

0 

V 

u 

IS 

41 

0 

9 

0 

0 

42 

4 

0 

2  J 

26 

0 

40 

i 

30 

43 

0 

7 

0 

0 

39 

7 

0 

24 

a 

0 

«0 

1 

•  TOTAL 

202 

28 

laa 

42 

.'14 

16 

I9« 

34 

U4 

i]i 

IK 

SI 

•  Total  of  the  5  durations  for  stations  in  this  region 

Our  The  duration  the  data  set  represents  le  average  JO  day  low  flow 

Sig  Iht  nuinbtr  uf  Mationt  whlth  Mmm   signiruant  Jtp«nutni«.  Iiind.  non  nnocaintss 

Not  The  nunoar  uf  tlitiuns  snun  \imm    wmepeniitnit .  ''••  'rua  trend.  ii<o  mtoiaBitit 

Per  in*  percent  binomial  proMDilKy  that  this  ntaotr  of  stations  oouid  ran  the  non  p«iarMl' 


Tables  3  and  4  show  the  results  of  data  from  >  20  years  of  record 
with  Table  3  showing  less  non-regulated  stations  than  the 
regulated  stations  on  Table  4.  The  regulated  stations  show  a 
greater  percentage  of  values  which  can  be  said  to  "pass"  the 
non-parametric  test  compared  to  those  which  "failed".  Although 
the  Central  Region  data  is  showing  computed  randomness  at  the  17. 
level  of  significance,  the  binomial  probability  is  still  "0"  for 
al 1  other  tests. 

Tables  5  and  6  summarize  the  analysis  results  for  stations  with 
less  than  20  years  of  record.  The  findings  suggest  that  neither 
the  length  of  record  or  the  effect  of  regulation  can  account  for 
the  overall  conclusion  that  the  statistical  tests  indicate  some 
degree  of  dependence,  trend  and  non-randomness  of  the  low  flow 
data  sets. 

It  is  possible  that  application  and  Interpretation  of  the 
non-parametric  tests  at  the  10%  level  of  significance  could 
result  in  fewer  statons  "failing"  the  non-parametric  tests. 
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Randomness 

Day 

IX 

5X 

IX 

5X 

IX 

51 

Duration 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

1 

6 

0 

0 

6 

0 

0 

5 

0 

2 

0 

2 

0 

6 

0 

0 

3 

6 

0 

0 

6 

0 

0 

5 

0 

3 

0 

2 

0 

6 

0 

0 

7 

6 

0 

0 

6 

0 

0 

5 

0 

3 

0 

2 

0 

6 

0 

0 

15 

6 

0 

0 

6 

0 

0 

5 

0 

2 

0 

2 

0 

6 

0 

0 

30 

6 

0 

0 

6 

0 

0 

5 

0 

2 

Q 

2 

0 

6 

0 

0 

•  TOTAL 

30 

0 

30 

0 

25 

5 

18 

12 

10 

20 



30 

0 

Central  Region 


Independence 

Trend 

Randomness 

Day 

U 

5X 

IX 

5X 

IX 

5X 

Duration 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

1 

3 

0 

0 

8 

0 

0 

8 

0 

0 

6 

2 

0 

0 

9 

0 

8 

0 

3 

a 

0 

0 

8 

0 

0 

7 

1 

0 

7 

1 

0 

0 

9 

0 

8 

0 

7 

8 

0 

0 

7 

1 

0 

7 

1 

0 

6 

2 

0 

0 

9 

0 

8 

0 

15 

8 

0 

0 

7 

1 

0 

7 

1 

0 

7 

1 

0 

0 

9 

0 

8 

0 

30 

7 

1 

0 

7 

1 

0 

7 

1 

0 

7 

1 

0 

0 

9 

0 

8 

0 

•TOTAL 

39 

1 

37 

3 

36 

4 

33 

7 

0 

45 

40 

5 

Southeastern  Region 


Independence 

Trend 

Randomness 

Day 

IX 

5X 

IX 

5X 

IX 

5X 

Duration 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

Sig. 

Not 

Per. 

1 

12 

0 

0 

12 

0 

0 

12 

0 

0 

0 

3 

0 

0 

3 

12 

0 

0 

12 

0 

0 

12 

0 

0 

0 

8 

0 

0 

7 

12 

0 

0 

12 

0 

0 

12 

0 

0 

0 

8 

0 

0 

15 

12 

0 

0 

12 

0 

0 

12 

0 

0 

0 

8 

0 

0 

30 

12 

0 

0 

12 

0 

0 

12 

0 

0 

0 

8 

0 

0 

•  TOTAL 

60 

0 

60 

0 

60 

0 

55 

5 

40 

25 

60 

5 

Dur 
Sig 
Not 
Per 


Total  of  the  5  durations  for  stations  in  this  region 

The  duration  the  data  set  represents  ie  average  30  day  low  flow 

The  nuneer  of  stations  which  show  significant  dependence,  trend,  non  randomness 

The  nuneer  of  stations  which  show  Independence,  free  from  trend,  and  randomness 

The  percent  Binomial  proDaoillty  that  this  nuneer  of  stations  would  fall  the  non  parametric  tests 


